esoph <- data.frame(
 agegp = structure(ordered(c(1, 1, 1, 1, 1, 1,
	1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
	2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4,
	4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5, 5,
	5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6), levels=1:6),
	class = c("ordered", "factor"),
	.Label = c("25-34", "35-44", "45-54", "55-64", "65-74", "75+")),
 alcgp = structure(ordered(c(1, 1, 1, 1, 2, 2, 2, 2,
	3, 3, 3, 4, 4, 4, 4, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4,
	4, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 1, 1, 1, 1,
	2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 1, 1, 1, 1, 2, 2, 2, 3, 3,
	3, 3, 4, 4, 4, 4, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4, 4), levels=1:4),
	class = c("ordered", "factor"),
	.Label = c("0-39g/day", "40-79", "80-119", "120+")),
 tobgp = structure(ordered(c(1,
	2, 3, 4, 1, 2, 3, 4, 1, 2, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3,
	4, 1, 2, 3, 4, 1, 2, 3, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1,
	2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2,
	3, 4, 1, 2, 3, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 4, 1, 2, 3, 4, 1,
	2, 1, 2), levels=1:4),
	class = c("ordered", "factor"),
	.Label = c("0-9g/day", "10-19", "20-29", "30+")),
 ncases = c(0, 0, 0, 0, 0, 0, 0, 0,
	0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 3, 1, 0, 0, 0, 0, 0, 2, 0,
	2, 1, 0, 0, 0, 6, 4, 5, 5, 3, 6, 1, 2, 4, 3, 2, 4, 2, 3, 3, 4,
	9, 6, 4, 3, 9, 8, 3, 4, 5, 6, 2, 5, 5, 4, 2, 0, 17, 3, 5, 6,
	4, 2, 1, 3, 1, 1, 1, 1, 2, 1, 2, 1, 0, 1, 1, 1, 2, 1),
 ncontrols = c(40, 10, 6, 5, 27, 7, 4, 7, 2, 1, 2, 1, 1, 1, 2, 60, 14,
	7, 8, 35, 23, 14, 8, 11, 6, 2, 1, 3, 3, 4, 46, 18, 10, 4, 38, 21, 15, 7,
	16, 14, 5, 4, 4, 4, 3, 4, 49, 22, 12, 6, 40, 21, 17, 6, 18, 15,
	6, 4, 10, 7, 3, 6, 48, 14, 7, 2, 34, 10, 9, 13, 12, 3, 1, 4,
	2, 1, 1, 18, 6, 3, 5, 3, 3, 1, 1, 1, 2, 1)
)
